Multiferroic materials for spin-based logic devices 
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Logical devices based on spin waves offer the potential to avoid dissipation mechanisms that limit 
devices based on either the charge or spin of mobile electrons. Multiferroic magnetoelectrics, which 
are materials that combine ferroelectric and magnetic order, allow direct switching of magnetic 
order and thence of spin-wave properties using an applied electric field. The intrinsic coupling 
between polarization and magnetic moments, generated by strong electronic correlations in these 
multiferroic materials, is argued to provide new approaches to spin-wave injection and spin-wave 
switching using applied voltages with no external magnetic field. These effects are shown to arise in 
a phenomenological Landau theory of coupled electronic and magnetic orders in multiferroic BiFeOa, 
and found to depend subtly on differences between the crystalline and film states of this material. 



I. INTRODUCTION 

Some limitations of current semiconductor logic de- 
vices, such as their relatively high power dissipation per 
operation, may be removed in devices that use electron 
spin rather than electron charge to carry out logic op- 
erations. While ordered domains of electron spin are 
currently the basic element of most information storage 
devices, information manipulation (i.e., logic operations) 
is currently accomplished using electron charge rather 
than spin. Spin-based logic is a relatively new field and, 
at the present time, has not convincingly demonstrated 
advantages in practice over ordinary charge-based logic. 
In many cases, it seems that devices using spin currents 
carried by itinerant electrons suffer from the same dissi- 
pation mechanisms as ordinary electronic devices. 

However, ferromagnetic and antiferromagnetic insula- 
tors provide another route to spin currents: they support 
"spin waves" that can transport spin and hence informa- 
tion over macroscopic distances, without accompanying 
transport of charge. This article discusses the relevant 
physics of "multiferroic magnetoelectric" materials that 
allow direct electrical control of the spin-wave propaga- 
tion, at room temperature and without external magnetic 
fields. Recently considerable effort has gone into demon- 
strating basic logic operations with spin wavesi^, as the 
dissipation mechanisms in this approach to logic and the 
interactions between spin waves are not well understood. 
Multiferroic materials with both magnetic and ferroelec- 
tric order address several of the most pressing problems in 
making such spin- wave devices practical for applications: 
in the near term, they allow electrical creation and de- 
tection of spin waves in a compact volume and can hence 
serve as efficient converters between electronic and spin- 
tronic logic. In the longer term, multiferroic materials 
may enable new designs for spintronic logic devices. 

The value of multiferroic magnetoelectrics is that these 
strongly correlated materials couple strongly to applied 
electric fields, because of their ferroelectric character, 
while at the same time showing magnetic order that en- 
ables propagation of spin waves. Furthermore, the mag- 



netic and ferroelectric orders are strongly coupled, as in 
the recent demonstration of room-temperature switch- 
ing of magnetization direction by an electrical field in-. 
As an example of how this enables spin-wave manipu- 
lation, we will show that the properties of a spin wave 
at wavevector k, such as its group velocity, depend on 
the angle between k and the Neel magnetization direc- 
tion. The result of this combination of effects is that elec- 
trical switching of the magnetization direction modifies 
the spin-wave propagation sufficiently strongly to stop or 
pass spin waves depending on the electrical state. This 
functionality is similar in concept to the function of fer- 
roelectric materials as electrooptic modulators-. While 
multiferroic materials have a relatively long history^i^, 
recent advances in fabrication and control of these mate- 
rials have led to a considerable resurgence of interest^. 

The theoretical calculations described in this paper 
are based primarily on a time-dependent Ginzburg- 
Landau model for the free energy of multiferroic ma- 
terials. The dynamical equations combine the familiar 
Landau-Lifshitz dynamics of ferromagnetic and antiferro- 
magnetic materials, the standard macroscopic free energy 
for polarization modes in ferroelectric materials, and ad- 
ditional terms that couple the two orders. The equations 
are solved to yield the spectra of coupled polarization 
and spin-wave modes. In choosing the microscopic de- 
tails of the multiferroic couplings, we focus for simplicity 
on the material BiFeOs. In bulk crystals of BiFeOa, cy- 
cloidal antiferromagnetic order is established at 650 K, 
well below the ferroelectric transition at 1120 K. The 
next section introduces the basic model for these crystals 
and explains one key result for spin-wave manipulation 
in these crystals well below the transition temperature: 
a zero-wavevector electric perturbation, such as an AC 
electric field, couples to finite- wavevector spin waves, me- 
diated by the cycloidal magnetic structure in the ground 
state. Films of BiFeOs differ significantly from bulk crys- 
tals; they show uniform antiferromagnetic order with a 
weak ferromagnetic component, rather than cycloidal or- 
der. Since the ability to create films is of considerable im- 
portance for applications, Section III gives an overview of 
how the model and results need to be modified for BiFeOa 



films, and concentrates on how the switching process al- 
ready demonstrated in one experiment'' will modify spin 
wave propagation. While our focus will be on BiFeOa 
because it is the predominant room-temperature multi- 
ferroic material, similar models can be applied to other 
low-tcmpcrature multiferroic materials as well. 



II. MAGNETOELECTRIC COUPLING IN BULK 

CRYSTALLINE BIFEO3: ELECTRONICALLY 

ACTIVE SPIN WAVES 



The first step in the theoretical description of this ma- 
terial is to construct an appropriate phenomenological 
free energy. At room temperature, the system is far be- 
low its critical temperature for antiferromagnetic and fer- 
roelectric transitions, so that thermal fluctuations can be 
ignored in the dynamics. The long-period incommensu- 
rate cycloid is produced in our model^ by a Lifshitz term 
in the effective Ginzburg-Landau free energy; the exis- 
tence of such a term was previously argue d^'^^i^^'^^ as the 
simplest symmetry-allowed explanation for the cycloidal 
order that is observed in experiments unless a very large 
magnetic field is applied (around 18 Teslapi^). The model 
free energy is 



orders are given by 
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Here L — \ Mi — Af 2 1 is a Neel vector describing the stag- 
gered sublattice magnetization, M = |Afi + 7Vf2| is the 
total magnetization of the material, and Pz is the mag- 
nitude of the ferroelectric polarization along z [the cubic 
(111) and equivalent directions in BiFeOs]. This free en- 
ergy is believed to capture the most relevant physics of 
BiFeOa bulk crystals, while for films, the dominant cou- 
pling between magnetic and ferroelectric orders is altered 
(see below). 

If the Lifshitz term had been omitted, the AFM and 
FE orders would be decoupled, leading to a simple ground 
state that is a commensurate, isotropic Heisenberg anti- 
ferromagnet with |Lop = ^, plus an easy- axis ferro- 
electric with polarization Pq = PqZ, Pq — ~a/u. Clearly 
A < and r > for an antiferromagnetic ground state, 
while a < and a± > for a ferroelectric ground state. 
The Lifshitz term induces a cycloidal ordering: the anti- 
ferromagnetic moment rotates in an arbitrary plane in- 
cluding z. It also increases the magnitude of Lq and 
Po . For a reference ground state about which to consider 
small perturbations^, we take a cycloid with 



Lo(a;) = Lq [cos{qx)z + sin(ga;)x] . 



(2) 



The cycloid direction is parallel to x and the pitch is 
q = aPo/c' . The strength of antiferromagnetic and polar 
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The linearized equations of motion for the ferroelectric- 
antiferromagnet order parameters include the familiar 
Landau-Lifshitz equations, complemented by the Debye 
equation of motion for the ferroelectric order parameter: 
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The constant 7 is a gyromagnetic ratio with dimensions 
of (sG)^^, while / has dimensions of s^^ and plays a 
similar role in the ferroelectric equation of motion. 

Details of the solution of these equations can be found 
elsewhere, along with a description of their coupling to 
external electromagnetic fields'^; in this work we would 
like to concentrate on the consequences for the prospec- 
tive applications discussed in the Introduction. The main 
result of interest for applications is that an oscillating 
electric field at zero wavevector can couple to spin waves 
at nonzero wavevector, as a consequence of the cycloidal 
order present in the ground state. At lowest order and 
with a perfectly harmonic cycloid, the AC electric sus- 
ceptibility has poles at the following shifted phonon and 
antiferromagnetic frequencies: 
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where contributions of order O {£,c'q^ /a±) were 
dropped. In these formulas, we have defined ^ ~ 
r(7Lo)V/ (~ 10 in BiFeOg) and the AC electric field 
is Ec^^'^*. The primed frequencies are defined in units of 
{jLo^/¥) - 10^2 rad/s. 

There are additional magneto-optical resonances that 
appear if a more realistic theory is constructed by includ- 
ing, for example, anharmonicity in the cycloidal ground 
state. The possible impact of these optically active spin 
waves for applications can be simply stated: the cy- 
cloidal structure of the material means that even zero- 
wavevector perturbations, which are potentially simpler 
to generate than the stripline-launched waves in cur- 
rent experiments^, can couple to nonzero- wavevector spin 
waves with a preferred direction. Since antiferromagnetic 
spin waves generically have large group velocity even at 
large wavelength, this may provide a viable means for 
efficient launching of spin waves with electrical fields. 

The cycloidal nature of the magnetic order leads to 
yet another potentially useful property for applications: 
The lowest frequency branch for spin wave propagation is 
strongly anisotropic. This occurs due to the presence of 
a non-trivial symmetry operation, namely we may rotate 



the cycloid pitch q out of the xz plane with no energy 
cost. This property implies the presence of a zero fre- 
quency mode with dispersion independent of ky-. 
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where we dropped contributions O(fc^). This soft mode 
dispersion is strongly anisotropic, and may be useful for 
electrical control of the spin wave group velocity via elec- 
trical switching of the Pq direction. Note that the spin 
wave group velocity is quite fast for propagation along 
the xz plane, which is here defined by the polarization 
Pq and the cycloid direction q. For this case the group 
velocity is given by ^LQ\frd ^ 10® cm/s. Interestingly, 
the anisotropic dispersion found in Eq. ([8]) leads to zero 
group velocity (at fc — > 0) for propagation along the y 
direction (perpendicular to the cycloid plane). 

The next section discusses the surprisingly different 
physics that occurs in BiFeOa films, where there is a uni- 
form canted antiferromagnetic ground state rather than 
a cycloid, in order to predict how electrical switching 
of the magnetization direction will modify the spin- wave 
propagation. 



III. SPIN- WAVE SWITCHING IN BIFEO3 
FILMS 

The cycloid order is destroyed in epitaxially constraned 
BiFeOa films, as evidenced by the measurement of a much 
larger magnetic susceptibilityi^ and confirmed by neu- 
tron scattering experiments^^. For this reason, the Lif- 
shitz invariant assumed in Eq. ([1]) can not be present or 
is too weak to produce noticeable effects (for example, 
the cycloid period in films is probably much larger than 
the typical domain size). Nevertheless, the films seem to 
have a weak macroscopic magnetization, of the order of 
0.02/ZB/Fe. This occurs even in the absence of applied 
fields and can only be explained by canting of the sub- 
lattice magnetizations. These considerations suggest the 
following model free energy for BiFeOs films: 
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Here Mi and M2 are sublattice magnetizations, and P is 
the ferroelectric polarization. We assumed^^, that cant- 
ing occurs due to a linear magnetoelectric effect through 
a Dzyaloshinskii-Moriya vector equal to dP. As a result, 
this free energy is invariant under inversion at a center lo- 
cated in between the two sublattices. This coupling is ex- 
pected if the low-temperature phase is imagined as result- 
ing from symmetry breaking in a high-temperature cubic 



phase, but another possibility has been suggested based 
on electronic structure calculationsiS. That work finds 
that the dominant contribution to the Dzyaloshinskii- 
Moriya vector is independent of P. In this case, the ferro- 
magnetic vector M is nearly decoupled from P. Exper- 
imental tests of how polarization switching modifies the 
magnetic structure in films are a high priority; tests so 
far indicate that the antiferromagnetic vector L switches 
to remain perpendicular to P, but the response of M 
has not been probed yet. 

The reference ground state for film structures is an 
uniform-canted state, 

Moi = Afo(sin/3* + cos/?y) , (10a) 

M02 = -Mo(-sin/?* + cos/3y), (10b) 

with canting angle approximately given by tan/3 ~ 
dPo/{2J). The spin waves for this canted state were cal- 
culated a long time ago for propagation along the sym- 
metry axis (Pq in our case), perpendicular to the cant- 
ing plane^^. There are two different spin wave modes: 
One is a low frequency mode, with linear dispersion lu = 
y/2Jc'k, where here a) — uj/{'^Mo). The other is a high 
frequency mode, with gap equal to the Dzyaloshinskii- 
Moriya coupling, J) w dPf) . 

We now discuss the possibility of electronic excitation 
of spin waves in a canted multiferroic such as BiFeOs 
films. The canted-uniform magnetic order of the films 
imply that an electrical probe with wavevector k will only 
excite spin waves at the same wavector k. This is in dras- 
tic contrast to the case of bulk BiFeOa described above, 
where the cycloidal magnetic order allowed the transduc- 
tion of a A: = excitation into an nq excitation. 

However, the presence of a macroscopic weak magne- 
tization in the film state allows a whole new set of mag- 
netoelectric functionalities. Below we provide a sum- 
mary of our results, with a focus on electrical control 
of magnon propagation. A detailed description is avail- 
able elsewhere^^. For example, the fact that the weak 
magnetization scales linearly with the electric polariza- 
tion (2Mosin/3 « dPoMo/J) allows electrical excitation 
of low frequency spin waves provided the AC electric field 
is polarized along the weak magnetization direction [x di- 
rection in our case). Moreover, an electric field polarized 
in the yz plane (perpendicular to the weak magnetiza- 
tion) will excite a high frequency magnon. The latter 
effect may be observed optically in the far infrared re- 
gion, since this will lead to an electromagnon resonance 
(a pole in the dielectric function) at the Dzyaloshinskii- 
Moriya frequency. 

Moreover, the presence of a macroscopic weak magne- 
tization will also lead to a sizable anisotropy in the low 
frequency magnon dispersion. The physical origin of this 
anisotropy is completely different than the one occurring 
in bulk BiFeOs [See Eq. ^ above]. In films, the pres- 
ence of a weak magnetic moment implies the spin waves 
will produce a macroscopic AC magnetization with finite 
amplitude SM. From Maxwell's equations we see that 



this must create a magnetostatic field of tlie form 



IV. CONCLUSIONS 
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where n is the spin wave propagation direction, k = kfi. 
The magnetostatic field contributes a term 2-k{5M ■ fif' 
to the free energy, tending to increase the spin wave fre- 
quencies whenever the quantity 5M has a finite projec- 
tion along n. As a consequence, the low frequency mode 
becomes anisotropic according to 
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Interestingly, a gap of the order of ■y/47r/J((iPo) appears 
for spin waves propagating along the y direction (the di- 
rection of the Neel vector in our case) . Similarly to the 
case of bulk BiFeOs, this effect may be used as an elec- 
trical switch for spin waves propagating along the xz di- 
rection. Changing the orientation of Pq electrically will 
drastically reduce the group velocity of a spin wavepacket 
propagating along the xz plane. 

The physical origin of this gap lies in the anticross- 
ing between the photon dispersion (uj = ck) and the 
low frequency magnon dispersion, a phenomenon anal- 
ogous to the formation of polaritons for optical phonon 
excitations. The magnetostatic propagation anisotropy 
Eq. (|12p arises precisely because the spin waves travel 
with finite k > lu/c, which is always greater than a few 
cm~^ due to the small domain sizes. For this reason, the 
magnon dispersion is effectively gapped, but we empha- 
size that the strict fc ^ limit, observable only at very 
long wavelengths (several centimeters) has w — > and no 
orientation dependence as expected by symmetry. 



In summary, we discussed the possibility of electrical 
excitation and control of spin waves in BiFeOs monocrys- 
tals and hcterostructures. The presence of the inhomo- 
geneous cycloidal order in bulk BiFeOs implies an AC 
electric field at fc « will excite spin waves at integer 
multiples of the cycloid wavevector (at k = nq). This 
shows that bulk BiFeOa is potentially useful as a source 
for finite wavevector magnons, which may become a vi- 
able alternative to the current strip-line designsi^. Al- 
though BiFe03 films do not have this functionality, spin 
waves may still be excited electrically, due to the presence 
of a Dzyaloshinskii-Moriya linear magnetoelectric effect. 

We remark that changing the direction of the ferro- 
electric polarization may force a switch of the Neel vec- 
tor orientation in a special class of multiferroic materials, 
an effect already demonstrated in BiFeOa films'^. In our 
work we showed in addition that BiFeOa monocrystals 
as well as hcterostructures possess a sizable anisotropy 
in the lowest frequency spin wave dispersion. The com- 
bination of these two physical phenomena implies that 
the magnon group velocity may be stopped electrically 
in BiFeOs, an important new functionality for devices 
based on spin wave propagation. Similar effects may be 
demonstrated in other well studied multiferroic materi- 
als, such as TbMnOg^, or BaNiF^^i. 
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